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We compare relativistic approximation methods, which describe gravitational instability in the 
expanding universe, in a spherically symmetric model. Linear perturbation theory, second-order 
. . . perturbation theory, relativistic Zel'dovich approximation, and relativistic post-Zel'dovich approxi- 

mation are considered and compared with the Lemaitre-Tolman-Bondi solution in order to examine 
the accuracy of these approximations. We consider some cases of inhomogeneous matter distribu- 
tion while the homogeneous top-hat model has been usually taken in the previous Newtonian works. 
It is found that the Zel'dovich-type approximations are generally more accurate than the conven- 
^ , tional perturbation theories in the weakly nonlinear regime. Applicable range of the Zel'dovich-type 

O ■ approximations is also discussed. 

^ , I. INTRODUCTION 

(N 

Structure formation in the universe is an important subject of research in cosmology. A standard view of the 
. structure formation is that density fluctuations with small ampUtudes in the early universe have grown to be a variety 
0^ of cosmic structures due to gravitational instability. The growth of the density fluctuations has been throughly 



investigated by linear perturbation theory of the Friedmann-Lemaitre-Robertson- Walker (FLRW) universe within 
both the Newtonian theory and general relativity 0|. Relativistic linear perturbation theory was flrst derived by 



I 

o ■ o 

\-[ \ Lifshitz Such relativistic treatments are indispensable when we consider large-scale fluctuations. In his theory, 

^ ' however, there remains a gauge problem that unphysical perturbations are included in the solutions. This problem 

^ ' was carefully studied later by Press and Vishniac Also developed was the gauge-invariant formulation Q, which 

•I— I : ^ 

. gives a conceptually straightforward way for dealing with cosmological perturbation. 
^ < It is true that the linear theories play an important role in the study of gravitational instability, but they are 
valid only in the region where density contrast 5 = {p — pb)l pb is much smaller than unity, [p is energy density of 
the perturbed FLRW universe and pf, is that of the background FLRW universe.) As 5 grows to be comparable to 
unity, nonlinear effects become essential and we need some kinds of nonlinear approximations. Tomita |^ developed 
second-order perturbation theory by extending Lifshitz 's work to study nonlinear effect of gravitational instability 
for the matter-dominated universe. His approach, however, still depends on the assumption of 5 being small. An 
approximation scheme without the assumption was proposed by Zel'dovich Q within the Newtonian framework. This 
scheme is known as Zel'dovich approximation, which is now widely applied to the problems of the large-scale structure 
formation. It has been shown that the Zel'dovich approximation can be regarded as a subclass of the first-order 
solutions in the Lagrangian perturbation theory Then higher-order extension of the Zel'dovich approximation, 
say, post-Zel'dovich approximation (and post-post-Zel'dovich approximation and so on), is straightforwardly derived 
via the higher-order Lagrangian approach |p|-pO|. Relativistic versions of the Zel'dovich approximation have been also 
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studied for the last few years by several authors [[ll|-|14[. Here we will focus on our tetrad-based approach, whose 
correspondence to the original Zcl'dovich approximation is made clear in Ref. JT2| and extension to second order is 
presented in Ref. |p^ . 

One of remarkably advantageous points of the Zel'dovich-type approximations, both the original Newtonian one and 
the relativistic version, is that they include exact solutions when the deviation from the background FLRW universe is 
locally one-dimensional. These exact solutions are known as Zel'dovich solutions ||^ in the Newtonian case and (some 
class of) Szekeres solutions [|l^,|6| in the general relativistic case, respectively. For this reason, the Zel'dovich-type 
approximations are presumably accurate in description of nearly one-dimensional collapse. It is not clear, however, 
whether they also give high accuracy in the case of non one-dimensional collapse. In the Newtonian framework, it has 
been investigated by using spherical models: The so-called top- hat collapse model the top- hat void model [p^ . 



and some more general case [|10[. (See also Ref. |19 for review.) In addition, there is also a recent work in which 
homogeneous spheroidal models are considered. An interesting implication is obtained in it: As the deviation of the 
models from the spherical symmetry becomes larger, the accuracy of the Zel'dovich-type approximations increases 
while the conventional (Eulcrian) approximations have the opposite tendency. It indicates that the Zel'dovich-type 
approximations may be the least accurate in the exactly spherical case. Then, considering the spherical case may tell 
us the lowest accuracy of the Zel'dovich-type approximations. 

In general relativity, an exact solution of the spherically symmetric dust model is known as the Lemaitre-Tolman- 
Bondi (LTB) solution It is, therefore, of interest to test the Zel'dovich-type approximations with the exact 

solutions to examine accuracy of the approximations. It is also instructive to clarify the differences between the 
conventional perturbation theories and the Zel'dovich-type approximations by imposing spherical symmetry and 
comparing them with the exact solutions. 

In this paper, we compare the relativistic approximations such as the linear perturbation theory by Lifshitz 
the second-order perturbation theory by Tomita Q, the relativistic Zel'dovich approximation by Kasai and 
the relativistic post- Zel'dovich approximation by Russ et al. with the LTB solution. We consider some initial 
conditions which have inhomogeneous matter distribution, not homogeneous one like the top-hat model adopted in 
the Newtonian works. It will be shown that the Zel'dovich-type approximations arc more useful than the conventional 
ones in the quasi-nonlinear regime in each case. 

The plan of this paper is as follows. In the next section, we summarize the relativistic perturbation theories 



mentioned above. In the section III. the LTB solution is introduced and the relation to the relativistic perturbation 
theories is considered. Main results of this paper are shown in the section IV . The section ^ contains summary of 
our results and discussions. 

Throughout this paper, units are chosen so that c = 1. Indices ^^v, - ■ ■ run from to 3 and ■ ■ ■ run from 1 to 3. 

II. RELATIVISTIC PERTURBATION THEORIES 

In this section, we summarize general relativistic perturbation theories which describe gravitational instability in 
the matter-dominated FLRW universe. We assume that the background is the Einstein-de Sitter spacetime, whose 
line element is 

ds^ = -dt^ + a^{t) {dR^ + R^de"^ + R^ sin^ Odct)'^) = ~dt^ + a^{t) kij dx'dx^ , (2.1) 
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where a{t) = t^/^ is the scale factor. The background four-velocity and energy density of the matter are = (1, 0, 0, 0) 
and pfc = \/{QiTGt^), respectively. Density contrast 5 shown later is defined as 6 = (p — pb)/ pi,, where p is energy 
density of the perturbed FLRW universe. 



A. Conventional linear and second-order theories 



Lifshitz pioneered linear perturbation of the FLRW universe in the synchronous gauge 

ds^ = —dt^ + Oil dx^dx^ . 



(2.2) 



In his theory, the scalar-mode solutions for pressureless matter (dust) are 



= 0, 



(1 



2t4*|. 



2t-i$ 



"(1) 
6(1) 



(2.3) 



Ifc " ' \k ' 

where ^' = ^'(a;) and $ ~ ^{x) are spatial arbitrary functions of first-order smallness, | denotes the covariant 
derivative associated with the background three- metric kij , and u' represents spatial component of the four- velocity 
of the matter. Subscript (1) denotes first-order perturbation quantity. We will not consider contribution of the 
decaying mode, which is proportional to t~^, later on. 

Tomita ||^ developed second-order perturbation theory by extending Lifshitz's work. He obtained the following 
second-order perturbative solutions from the first-order scalar mode solutions: 



7y- 
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'5(1) + ^(2 



0, ^ 
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0, 
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(2.4) 



Subscript (2) represents second-order perturbation. Here we neglected the second-order tensor mode, which is induced 
by the first-order scalar mode and does not appear in the spherical case. 



B. Zel'dovich-type approximations in general relativity 

In this subsection, we review a relativistic version of the Zel'dovich approximation developed by us [p^|jl3| . The 
irrotational dust model is assumed and then we can take the comoving synchronous coordinate 

ds^ ^ -dt^ + gi-j dx'dx^ (2.5) 

with the four-velocity = (1,0,0,0). Thanks to this choice of the gauge, the energy equation w^T^^^ = with 
rptj.f ^ diag[/9, 0, 0, 0], which becomes 

p + pK\ = 0, (2.6) 
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is formally solved in the form 



det [gij{Un,x) 



(2.7) 



det [gij {t, x)] 

Here an overdot (') denotes d/dt, and K^^ is the extrinsic curvature, whose expression in the present gauge is 

J 2i 



K^j — ^g^'^g'jk- Introducing the triad 



gij = a [t)6(^km eyey 



(2.8) 



Eq. (2.7) is rewritten as 



P = Pb 



det [e'i'{tin,x)\ 
det[e^f (i,a;)] 



(2.9) 



We obtain perturbative solutions for the triad e\ regardless of the energy density p up to second order in the 
following form 



(2.10) 



where k'^f is the background triad defined by fey — (5(fc)(^) fc'^J^''fc^^\ and E'^f and e^^"* are the first-order and the 
second-order solutions given by 



e(f = fc(f (ti + it ^\ 



(2.11) 



Here — '^{x) is the same function as the one used in Eqs. (2.3) and (2.4), and ip^ ^ = tp'^j{x) and (f^^ — (p^ ^{x) are 
quadratic quantities of written by 



(2.12) 



14 



(2.13) 



Note that we removed a remaining gauge freedom in the linear level to derive the above solution. (See Appendix A 
of Ref. (l^.) And we again neglected contributions of the decaying scalar mode and the tensor mode. We find that 
the solution ( ^.10 ) is consistent in the metric level with Eqs. (2^) and {2A) |p^ . 



Relativistic Zel'dovich and post-Zel'dovich approximations are obtained by substituting Eq. ( p.lC| ) into Eq. ( p. 9]) 

-1 



det 



-1, 



JPZA 



= det 



S' ■ H ^ — 

1 + f* 



(2.14) 



(2.15) 



Abbreviations ZA and PZA denote the Zel'dovich and the post-Zel'dovich approximations. 
As was written previously, the results of the conventional linear and second-order theories are 
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W = -i^*'lfc , (2.16) 



^SEC 



ti^\ + Iti (^I'^vi;,, + 6*vi/iy + (5(vi/if,)2 + 2*i;=,*y . (2.1?) 



Abbreviations LIN and SEC denote the linear and the second-order perturbation theories. Expanding Eqs. (2.14) 



and (2.15) under the condition \\^\\ <C 1 (where ||^'|| denotes an appropriate norm of a function '5), Eqs. (2.16) and 



( ^.17 ) can be also obtained, respectively. In this sense, ZA and PZA are extensions of LIN and SEC to \5\ ^ 1. 



III. SPHERICALLY SYMMETRIC MODEL 

In this section, we consider the spherically symmetric model of gravitational instability in the FLRW universe. 
There exists an exact solution known as the LTB solution, which includes three arbitrary functions, in the spherically 
symmetric case [^l| . Here we make clear the relations between the arbitrary functions included in the LTB solution 
and the ones which appears in the approximation methods mentioned in the section 0. The line element of the LTB 
solution is 

ds^ = -dt^ + dR" + r2(d6|2 + s\n^ 9 dcp'^) , (3.1) 

where ( ') = d/dR and / = f{R) is an arbitrary hmction which is related to initial velocity of dust, r = r{t,R) 
satisfies the following differential equation 

r' = ^ + f{R) (3.2) 



with an arbitrary function F{R), which represents initial distribution of matter. Eq. (3.2) can be integrated as follows 
(i) / > 0: 

F F 

r = —{coshrj-l) , t-to{R) = jjjj^{sinhr]-7]) , (3.3) 



(ii) / < 0: 



(iii) / = 0: 



F F 

:(l-cos77), t-to(i?)= (?7-sin?7), (3.4) 



-2f' ' 2(-/)3/2" 



9F\ ^ 2 

— j {t~to{R))^, (3.5) 

where to(^) is an integration constant. The above cases (i), (ii) and (iii) may be called "open," "closed," and "flat," 
respectively, as in the FLRW universe. In these three cases, the density reads 

8^Gp^^. (3.6) 

Apparently the LTB solution includes the three arbitrary functions, f{R), F{R), and tQ{R). But dynamical degree 
of freedom is actually two because there remains freedom of choice of the radial coordinate R. 

The LTB solution is some extension of the FLRW solution and is often used as a model of an inhomogeneous 
universe (see Ref. ||l^ for review), and this solution can be reduced to the FLRW solution. By choosing f — 
and to = (and F — ^R^ for convenience), we easily find that r = a(t)R and the LTB solution is reduced to 
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the spatially flat FLRW solution (2.1). Furthermore, we can represent this solution by the form of the flat FLRW 
solution with its perturbations and can see that the arbitrary functions in the LTB solution corresponds to those 
in the linear perturbation theory. To see this, we consider spherical linear perturbation of the spatially flat FLRW 
solution. Substituting 
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/ = /(!), ^0=^0(1), F = -i?-^ + , r = aR + r^i 



(3.7) 



into Eq. (3.2), the linearized equation for r(j^ can be obtained, where the quantities with subscript (1) are treated as 
linear perturbation. Solving this equation, we obtain 



r(i) = aR [ -F(i)R- 



20 



t-^B 



(3.8) 



where B = B{R) is an integration constant. Actually B is related to io(i) by i? = —"^tQi^iy It is easily seen by 



choosing f — and using Eq. (3.5), which reads 

^(1) 



aR { ^F(^^R-^ 



0(1) 



(3.9) 



Then we can write "linearized LTB metric" in terms of /(ij, <o(i)7 ^J^d F(j) as follows 

IBM = 1 + \{F(i)R-^y - /(I) + ^t^{hi)R-'y - ^i-^(<o(i)i?)' 



lee = i?' 1 



2^(1)^ ^■ 



10 



0(1) 



(3.10) 



On the other hand, the solution of the linear theory (2.S) gives 

20 2 



9 

R' (l 



20 



Comparing Eqs. ( 3.1C ) and ( 3.11 ), we find the following relations 

i^(l) = ^*i^^ = t^^,)^-^-^'R-' 



(3.11) 



(3.12) 



This tells us that the arbitrary functions f{R) and to(i?) correspond to the growing and the decaying modes, respec- 
tively, in the linear level. We choose to(-R) = hereafter because contributions of the decaying mode are not taken 
into account throughout. Moreover, we know that and /(i) are related to each other by the function . This is 
because we eliminated a residual gauge freedom and thus fixed the gauge condition completely in the previous section. 
This fixing corresponds to determination of the choice of the radial coordinate R in the spherical case. 

To see the relation between the LTB solution and the second-order perturbation, we introduce /(2) i F(2) i and r(2) 
so that 



20 , 4 „ 

/-y*'i? + /(2), ^=9^' 



40 
27 



^R^ + F(2) , r = aR + r(i) + r(2) 



(3.13) 



and make calculations in the same way. Here /(2), ^"(2); and r(-2) should be regarded as 0(||^|p) and r(2) is obtained 
by solving Eq. (3.2) perturbatively. Then we obtain the part of 0(||5'|p) in jrr and jeg. Comparing these jrr and 



7ee with the solution of the second-order theory (2^), we find 



6 



^(2) = ^*'i?', /(2) = ^(*"i?'-2**'i?). (3.14) 



On the other hand, the LTB solution as an exact model is obtained by choosing 



F = + + l^^vI/2i?3 (3 15) 

9 27 243 ' ^ ' 

f = ^^'i? + - 2**'i?) . (3.16) 



From Eq. (3^), the density contrast of the LTB solution reads 

3 F' 

^^TB = 4^-1. (3.17) 
Let us turn our attention to the approximation methods. If we impose — ^{R), Eqs. ( 2.16 ), ( ^.17 ), ( 2.14 ) and 



( P.15[) become 



5lin = -t^ (*" + 2*'i?-i) , (3.18) 
5sEc = -i* + 2*'i?-i) + ^ti + 6*(*" + 2*'i?-i)) + (^"^ + + ^*'2i?"2^ , (3.19) 



X 



1 + 2 ^ i7 ^ i\ 3 

1 + f* y ^ 



Peculiar velocity, which represents deviation of motion of dust shell from the Hubble expansion and is defined by 
V = f — Hr, (where H = a/ a is the Hubble parameter) is written as follows 

2 1 

VLIN = VZA = ^t^-^' . (3.22) 

VSEC--VPZA = hH' + ^tH-i>'^R-4^'^')-^m'^R-' . (3.23) 

Note that vlin = vza and vsec — vpzA because, in the metric level, the Zel'dovich-type approximations coincide 
with the conventional ones. 

Now the density contrast and the peculiar velocity of the LTB solution and the approximations are written in 
terms of the only one function The function '5 should be determined from initial conditions so that the regularity 
conditions at R = 0, i.e., ^'(i? = 0) = and ^''(i? = 0) = are satisfied. Then the peculiar velocity at i? = is 
always zero both in the LTB solution and the approximation. Moreover, if 5* is taken so that cx i?^ near R = 0, 
the peculiar velocity near i? = is proportional to i?, w cx R. 

IV. COMPARISON OF LTB SOLUTION AND APPROXIMATIONS 

Let us proceed to comparison of the LTB solution and the approximations. As mentioned in the section |l|, Sza and 
SpzA include Slin and 5 sec, respectively, when the density contrast is small. As for the peculiar velocity, ZA and 
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PZA are coincident with LIN and SEC, respectively. Moreover, expanding the density contrast of ZA in the following 
form 

10 



5zA ^ -t^i'i'" + 2^'R- 



9 



3^/2^-2 



)+0(||*in, (4.1) 



it is found that ZA includes second-order (and higher) terms partially in the expression of the density contrast. Thus 
we can expect that Sza is as accurate as Ssec late time. 

In order to investigate relative accuracy of the approximations quantitatively, we compare the LTB solution and 
the approximations by using some specific initial conditions. Here initial conditions can be completely fixed by giving 
an initial density profile din{R)- For simplicity, we assume that din{R) is a first-order quantity. Then the arbitrary 
function 5* is determined by the relation 



JLIN 



= -(*" + 2*'i?-i) = 4„(i?) 



(4.2) 



with normalization tin = 1- (tin may be regarded as the decoupling time in the history of the expanding universe.) 
Eq. (4.2) is solved to determine the function with the boundary conditions '!'(/?, = 0) = and ^''(i? = 0) = 0. 
Here we consider the following two cases: 

S^niR) = 

and 



R_ 



exp 



R 

Rq 



(4.3) 



Sin{R) = e 



R 

Ro 



1 + ^- ^ 



R 
Ro 



R 

exp 1 - — 



(4.4) 



where e is a small constant (|e| <C 1) which represents amplitude of an initial density perturbation, and Rq is a 
comoving scale of the fluctuations. The former case can be regarded as smoothing out the top-hat model while in the 
latter case, if e < 0, the neighborhood of i? = is underdense region (void) and the outside is overdense, and then 
the shell-crossing will occur. (If e > 0, the latter case also shows only similar behavior to the top-hat model as the 
former one.) For both of the cases, we can evaluate 6 at the center of the fluctuations {R — 0) from the LTB solution 
and the approximation methods in the following form 



^LTB 



(i? = 0) 



9 (?7 — sin 7])^ 

^ (1 - cosri)'-^ 

9 (?7 — sinhr])^ 

^ (cosh 77 - 1)^ 



with 



t = 




/ Q 3 

- (— e)^2 (sinh?7 — 77) 



for e > , 

for e < , 

for e > , 

for e < , 



(4.5) 



and 



(i? = 0) = et3 
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SsEc{R = 0) = et-^ +^eH-^ 



6za{R^0)= ( 1- 



-3 



JPZA 



(i? = 0) 



- 1 
^2 



■i3 



21 



-3 



(4.6) 
(4.7) 

(4.8) 
(4.9) 
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These expressions do not depend on the details of initial density profiles. It is of essence in the calculation that 
Sin — e near R — 0, and e > and e < correspond to / < and / > at i? = 0, respectively. The above results 



are described in Figures 1 and 2, where Eqs. (15)-(19) are plotted as functions of 5ltb{R — 0). Figure 1, which 
represents the collapse case, tells us that ZA is more accurate than LIN and PZA is more accurate than SEC in all 
5ltb{R = 0) > region. (Figure 1 shows only the region of < 5ltb{R = 0) < 1, but the tendency shown in Figure 
1 does not change in denser region 5ltb{R = 0) > 1.) We also find that ZA becomes more accurate than SEC at 
late time when 5ltb is larger than about 0.5. This is understood by comparing Eq. ( [4.7[ ) and the expanded form of 
Eq. K. 



5zA{R = Q)^eti + '^^t-^ +0{e'') . (4.10) 

(This form is also obtained from Eq. ([l.l|).) We see from Eqs. (4.7) and ( 4.10| ) that 5za is smaller than 6sec at early 
time [et^/'^ <C 1) due to the lack of the terms in O(e^). In this sense, 5za is less accurate than 5sec when et^/"^ <C 1. 
However, due to the existence of the singularity at et^/^ = 3 in 5za^ ^za becomes to be more accurate than Ssec a-t 
late time. This existence of the singularity in Sza essentially determines the asymptotic behavior of dzA- Indeed, the 



exact solution Sj^tb for e > in Eq. (4.5) has a pole of order three at et^^^ — 3(37r/2)^/^/5 ~ 1.7 (r/ — 2tt). This pole 
corresponds to the crunching time at i? = 0, and the singularity occurs at this time. 

From Eqs. ( |4.5| )-( |4j| ), we can also evaluate accuracy of the approximations quantitatively at the turnaround time 
(et2/3 = 3(97rV2)^/^/10 ^ 1.1), though it is not drawn in Fi gurc 1. Here the turnaround time is characterized 
by f = 0, i.e., the maximum expansion. Physically speaking, the density fluctuation begins to collapse due to 
gravitational instability, overcoming the cosmic expansion at the turnaround time. At the turnaround time, Sltb 
becomes to be 4.6. To this Sltb, Slin, Ssec, ^za and Spza grow to about 23%, 43%, 60% and 84%, respectively. 
It is natural that SpzA is more accurate than Sza from the viewpoint of the singularity at the crunching time. SpzA 
also has a pole of order 3 at et^/"^ = (\/133 — 7)/2 ^ 2.3. This crunching time is nearer to the real crunching time 

1.7 than that of Sza- 

From Figure 2, which denotes the void case, we find that PZA gives the best fit at early time before Sltb — —0.7 
while ZA works best at late time when Sltb is smaller than about —0.7. At late time, PZA gives bad results. It is due 
to difference of the signature between the first-order (-et^^'^/S > 0) and the second-order (-e^t^/^/21 < 0) terms when 
e < 0. The same feature also appears in SEC but at earlier time. However, in PZA, this difference is more serious 
than in SEC. In SEC, Ssec grows as < 0, while in PZA, the difference of the signature makes SpzA diverge 

at a finite time et^/"^ = (7-f \/l33)/2. This divergence is an apparent one which is caused by the formalization of PZA. 



Indeed, we can easily see from Eq. (4.5) that the exact solution Sltb has no pole in i > (i.e. there is no singularity) 
and approaches to —1 as ~ t^^. Though there is no singular point in Slin except t — cxd, Slin takes value smaller 
and smaller without limit as the time increases because there is no physics to stop this decrease of the energy density 
in this order. Then, only Sza predicts true asymptotic value of S without growth and apparent singularities. But Sza 
approaches to —1 as t^^. This difference is also seen in Figure 2. In the Newtonian case, detailed discussions on the 
void can be seen in Ref. [T^ . Our relativistic results up to now are quite similar to Newtonian ones, which are given 
in Refs. 0,|l^,|l|. 

Figures 3 (a) and 4 (a) give the density contrast as a function of R/Ro when e = 1.0 x 10~^ and t = 2.0 x 10^, 
and e = —1.0 x 10^"^ and t = 3.7 x 10** with the initial density profile (^^), respectively. These two figures show that 
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difference between the LTB solution and the approximations is the largest at i? = in the former case. Hence, it 
is sufficient to consider the difference at i? = when we examine the accuracy of the approximations in the former 
case (4.3). 



We also see the evolution of the pecuhar velocity with the initial density profile (4.3) in Figures 3 (b) and 4 (b). 
Here we consider the peculiar velocity normalized by the Hubble flow Hr. It will be convenient to use it to see the 
deviation of the model from the FLRW universe in the metric level. For example, the normalized peculiar velocity 
v/Hr — —1 at the turnaround time. Figures 3 (b) and 4 (b) show the normalized peculiar velocity corresponding 
to Figures 3 (a) and 4 (a). Although these figures show that the normalized peculiar velocity v/Hr is not zero at 
i? = 0, we must note that this is due to our normalization. As mentioned in the last section, the peculiar velocity 
must behave ~ R near R — Q. On the other hand, Hr also behaves R near R — Q. Hence our normalized peculiar 
velocity does not vanish at i? = due to the normalization. It is also noted that the peculiar velocity obtained 
from the Zel'dovich-type approximations is the same as the one obtained from the conventional approximations as 
mentioned in the previous section. Furthermore, we find from Figures 3 (b) and 4 (b) that the deviation from the 



FLRW universe is maximum near i? = for the initial profile (4.3). And this is also shows that it is sufficient to 



consider the difference at i? = when we examine the accuracy of the approximations in the former case (4.3) as 
mentioned above. 



On the other hand, for the initial density profile given by Eq. (4.4), it cannot be said that the largest deviation of 
the approximations from the LTB solution occur at i? = 0. Figures 5 (a) and 6 (a) are for this initial profile when 
e = —1.0 X 10^"^ and i = 2.0 x 10^, and e = —1.0 x lO^'' and i = 3.0 x 10^, respectively. For this initial density proffle, 
the shell-crossing singularity will occur. The tendency of the occurring of shell-crossing can be seen from the peculiar 
velocity in Figures 5 (b) and 6 (b), where the profile of the normalized peculiar velocities for this case are drawn. 
In Figures 5 (b) and 6 (b), there exists a u = point at R/Rq ^ 2.5. This point, which is denoted by R — Rc 
hereafter, is a boundary where the universe is locally "open" (/ > 0) and "closed" (/ < 0), and then / = at the 
point. The peculiar velocity in the void region R < Rc is positive, while that in the closed region R> Rc is negative. 
Then we can expect that the shell-crossing of the dust matter will form at R ^ Rc within a finite time. 

Indeed, one can see from Eqs. (|3.2|) and (p]^) that the shell-crossing, which is characterized by a finite radius at 



which r' vanishes will occur at the radius R = Rc- From Eq. ( |3.5D , which is the solution when / = 0, we know 

. (4.11) 



9^; 
4 



(Subscript c denotes value at i? = Rc.) Differentiating Eq. ( p.2[ ) with respect to R and using Eq. (4.11), we obtain 
the equation for r^. Integrating this equation, one finds 

'^'''''*»4(.^)*/;«'+Cri, (4,2) 



4 \9FcJ 5 \9Fc, 

where C is an integration constant and is not essential in our argument. Here, r' must be positive initially if r is a 
monotonically increasing function of R which is in our case. This means that the first term (cooperate with the third 



term) in Eq. (4.12) must dominate on the initial surface. However, since < at i? = Rc, r'^ must vanish within a 
finite time. Hence, at R = Rc, the shell-crossing singularity will occur. (More generic arguments about the occurrence 
of shell-crossing singularity can be seen in Ref. p3] and the shell-crossing may occur in the region / < at first.) 
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It should be noted that, in Figures 5 (a) and 6 (a), 5za and 5pzA take the same value as that of the LTB solution 
ai R — Rc where the shell-crossing will occur. Then we must say that these figures show that the Zcl'dovich- 
type approximations are not necessarily inaccurate even when the shell-crossing is occurring. Let us consider the 
reason here. In fact, the deviation from the background Hubble expansion is locally one-dimensional at the point. 
The definition of "locally one-dimensional deviation" we adopt here is that two of the eigenvalues of the peculiar 

H5 



deformation tensor 



J = ■^J^'^jjk are zero |23 . According to the definition, let us show the local 



one-dimensionality a.t R — Rc ^ 0- From Eqs. (pjj) and (4.11) 



^0, v\ 



R=Rc 



R=Rc 



R=Ro 



a 
a 



vltb 



(4.13) 



R.=Rc 



Since V^^ is diagonal in the spherical case, this means that the deviation at R — Rc is locally one-dimensional. It is 



known that the Zel'dovich-type approximations become exact when the deviation is locally one-dimcnsional |12| , p3[ . 
In this argument, the origin i? = must be excluded because all components of the peculiar velocity vanish at the 
origin. Thus we can lead a significant consequence that, at the points where / = 0, the Zel'dovich-type approximations 
coincide with the exact LTB solution, i.e., 

Sltb = SzA = SpzA and vltb = vzA = vpzA = . (4-14) 

Note that this consequence is not limited to the specific initial density profile (|4.3| ) nor ( [4.4| ). 

Turning to Figures 5 (a) and 6 (a), we see that the coincidence of the density contrast at R/Rq ^ 2.5 {R — Rc) 
contribute to accuracy of the Zel'dovich-type approximations, and 6za and SpzA give good fit around R/Rq ~ 2.5. 
This is the reason the Zel'dovich-type approximations do not necessarily give bad results even when the shell-crossing 
is occurring. 



V. SUMMARY AND DISCUSSIONS 



We have tested the relativistic perturbative approximations to gravitational instability with the LTB solution. 
It has been shown that the Zel'dovich-type approximations give higher accuracy than the conventional ones in the 
quasi-nonlinear regime |5| ~ 1 within general relativistic framework. Our results are partly similar to the Newtonian 
ones, but our consideration is more generic. Especially we considered some cases in which matter distribution is 
inhomogeneous, and found that the Zel'dovich-type approximations are not necessarily inaccurate even when the 
shell-crossing is occurring. Of course, the occurrence of the shell-crossing shows the break down of our treatment. 
However, this is due to the failure of our description of the matter as dust rather than the failure of Zel'dovich-type 
approximations. 

Indeed, one of the case considered in the previous section includes the f — point, where the universe is locally 
"open" inside and is locally "closed" outside and the shell-crossing will occur at this radius. We have seen, in general, 
at the f — points (except the origin R — 0), the deviation from the FLRW model is locally one-dimensional and 
the Zel'dovich-type approximations become exact. And in the neighborhood of the points, we can expect that the 
Zel'dovich-type approximations are particularly accurate. The case considered here is exactly such an example. It 
should be also noted that "one-dimensionality" which makes the Zel'dovich-type approximations exact means not only 
globally plane-symmetric but also locally one-dimensional: Such situations appear even in the spherically symmetric 
model. 
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To discuss applicable range of the Zel'dovicti-type approximations, we reconsider the density contrast at i? = 
in the collapse case. At the turnaround time, accuracy of the Zel'dovich-type approximations already begins to fall 
down, i.e., SpzA is about 84% of Sltb and 5za is about 60%. Inaccuracy will be accelerated beyond the turnaround 
time. In this sense, the turnaround epoch, when the peculiar velocity is as large as the Hubble expansion, is one of 
criterion of applicable range of the Zel'dovich-type approximations. However, this criterion might not be practical, 
because one cannot know the correct turnaround time in general, while we have been able to know that from the 
exact solution in our case. Instead of the turnaround time, the Zel'dovich-type approximations tell us the crunching 
time as the singularity in the density contrast, approximately. Furthermore, we have also seen that PZA tells us this 
crunching time more accurately than ZA in our spherical model. Then we may be able to know the approximate 
turnaround time by the half of this crunching time of PZA. 

It is said that the Zel'dovich approximation predicts pancake formation in the gravitational collapse of dust |^ . But 
it is beyond the turnaround time that the pancake will be formed and thus accuracy of the Zel'dovich approximation 
is not ensured at that time. If we try to examine final stage of the collapse quantitatively, we will need to develop a 
new approximation scheme which gives an accurate description even beyond the turnaround epoch. 
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Fig. 3 (b). Profile of normaUzed peculiar velocity corresponding to Figure 3 (a). 



16 




17 




1 2 3 4 5 



R/Ro 

Fig. 4 (b). Profile of normalized peculiar velocity corresponding to Figure 4 (a). 
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Fig. 5 (b). Profile of normalized peculiar velocity corresponding to Figure 5 (a). 
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